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Abstract. A search for neutrino-induced muons in correlation with a selection of 40 gamma- 
ray bursts that occurred in 2007 has been performed with the ANTARES neutrino telescope. 
During that period, the detector consisted of 5 detection lines. The ANTARES neutrino 
telescope is sensitive to TeV-PeV neutrinos that are predicted from gamma-ray bursts. No 
events were found in correlation with the prompt photon emission of the gamma-ray bursts 
and upper limits have been placed on the flux and fluence of neutrinos for different models. 
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1 Introduction 



Gamma-ray bursts (GRBs), transient flashes of gamma-rays having a duration of sub-seconds 
up to several hundred seconds, are the most powerful known extra-galactic events. Several 
models predict a burst of high-energy neutrinos in concurrence with the flash of gamma-rays, 
also referred to as the prompt emission (see e.g. [1] for a review). In this, neutrinos are 
produced by interactions of protons that are accelerated by shock waves with energetic am- 
bient photons. A variety of models has been put forth covering a wide range in the expected 
amount of neutrinos. Some of those models have already been challenged by the recent 
IceCube measurements [2]. In particular, the IceCube results do not support the notion of 
GRBs as primary source of ultra high-energy cosmic rays, at least under the assumption that 
these ultra high-energy cosmic rays are protons coming from the decay of photohadronically 
produced neutrons (i.e. via the resonant channel p + 7 — > A + — > n + tt + ) and are therefore 
connected to neutrinos. In this scenario, neutrons, as opposed to the protons, can escape 
because they are not magnetically confined inside the source. Nevertheless, a viable phase 
space of models with a potentially measurable neutrino flux from GRBs in the TeV-PeV 
region can still be tested for a better understanding of the nature of these energetic events. 

In this paper, the first data taken with the ANTARES neutrino telescope in 2007 are 
used to perform a search for TeV-PeV neutrinos in correlation with a selection of GRBs 
detected by satellite observatories. The data are treated in a stacking approach in which 
the data observed during the prompt emission of all the selected GRBs are accumulated. 
The detection of a single neutrino event would constitute an observation with more than 3a 
significance. 

Previous searches for neutrinos from GRBs have been performed in the Northern hemi- 
sphere by AMANDA [3] and IceCube [4, 5] at similar energies (TeV-PeV) and by ANITA [6] 
at higher energies (>PeV). IceCube also included some GRBs in the Southern hemisphere [2] 
that is observed by ANTARES but not during the period considered in this analysis. Other 
searches in the Southern hemisphere were performed by Super-Kamiokande [7] at lower en- 
ergies (MeV-100 TeV) and by the Baikal neutrino telescope NT200 [8]. Both these exper- 
iments have a lower sensitivity than the current analysis in the comparable energy range. 
The ANTARES neutrino telescope is the most sensitive instrument to observe high-energy 
neutrinos from the GRBs considered in this analysis. In addition, the vast majority of these 
GRBs have not been studied for neutrino emission in the TeV-PeV range before. Despite 
the fact that more stringent limits have been published [2, 4, 5] for different GRB samples, 
ANTARES might still have observed events from the GRB sample presented here. 

The completed ANTARES detector is the largest neutrino telescope on the Northern 
hemisphere and is sensitive to neutrinos in the TeV to PeV energy range. Located in the 
Mediterranean Sea, the ANTARES detector is sensitive to GRBs in the Southern hemisphere 
where the sensitivity of IceCube significantly suffers from the large background from muons 
produced by cosmic ray interactions in the atmosphere above the detector. Due to the tran- 
sient nature of the GRBs and the variety of their characteristics, it is essential to permanently 
monitor the full sky in order to maximise the probability to observe a neutrino signal. The 
search described in this paper is the first in a series of searches which will be pursued with 
the data that will be acquired by the completed ANTARES detector up to 2016 and then 
continued with data from the future KM3NeT [9] neutrino observatory in the Mediterranean 
Sea. 
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Figure 1. Schematic view of the ANTARES detector in its Enal configuration. 



2 Neutrino detection 

High-energy neutrinos can be detected indirectly by a neutrino telescope such as ANTARES. 
The detection principle relies on measuring the Cherenkov light induced by high-energy 
charged particles that are produced in a neutrino interaction inside or near the instrumented 
volume. In particular, a high-energy muon is produced in a charged current muon-neutrino 
interaction. High-energy muons can travel large distances, which facilitates an accurate 
determination of the direction. At the typical neutrino energies considered in this analysis 
(5 TeV-6 PeV), the direction of the muon closely follows that of the incident neutrino. Above 
5 TeV the median angle between the neutrino direction and the muon direction is less than 
0.3°. The direction of a high-energy muon can thus be correlated to the positions of GRBs. 

The Cherenkov light is detected by photo-multiplier tubes, housed in optical mod- 
ules [10]. Triplets of optical modules are attached to vertical detection lines of about 450 
metres height, which are anchored to the sea bed at a depth of 2475 metres and held upright 
by a buoy. In its final configuration, the ANTARES detector consists of 12 detection lines 
with a spacing of about 60 metres (see figure 1). Eleven of these detection lines each contain 
25 evenly spaced triplets of optical modules and one contains 20 triplets. The positions and 
orientations of the optical modules vary due to the sea currents. An acoustic positioning 
system, combined with compasses and tiltmeters located along the detection lines, measures 
the positions and orientations of the optical modules with an accuracy of about 10 cm. A 
detailed description of the detector is given in ref. [11]. 

The arrival time and charge of the photo-multiplier tube signals are digitised [12] and 
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transmitted to shore. The absolute time stamping is performed by interfacing the clock sys- 
tem to the GPS which provides a time accuracy of about 100 ns with respect to Universal 
Time Coordinated [13]. The data are dominated by optical background due to biolumines- 
cence and natural radioactive decays. On shore, the physics signals are selected from the 
data stream by a software data filter that operates in real time. This data filter has multiple 
algorithms implemented, each designed to find a specific physics signal. The physics events 
are stored on disk for offline reconstruction. A detailed description of the data acquisition 
and data filtering is given in ref. [14]. 

3 Data selection 

The analysis described in this paper has been applied to data collected in 2007, in the period 
from January 27 to December 7. At this time, the ANTARES detector was still under 
construction and consisted of 5 detection lines. 

3.1 GRB selection 

The GRBs examined for neutrino emission were selected from the observations of GRBs by 
satellite instruments, as archived by the Gamma-ray Burst Coordinates Network (GCN) [15, 
16]. In order to suppress the background from muons produced by cosmic ray interactions 
in the atmosphere above the detector, hereafter referred to as atmospheric muons, GRBs 
were selected that occurred below the horizon of the ANTARES detector. Neutrinos from 
such GRBs would traverse the Earth and cross the detector in an upward-going direction. In 
the period considered, 46 GRBs occurred below the ANTARES horizon during physics data 
taking. 

By requiring the availability of detector alignment data, and by applying quality criteria 
based on the environmental conditions, the data in coincidence with 6 GRBs were excluded. 
The time and position information of the 40 remaining GRBs, listed in table 1, were used 
to search for a correlated neutrino signal. A total of 32 GRBs were detected by the Swift 
satellite [17], 4 by INTEGRAL [18] and 4 by other satellites of the Third Interplanetary 
Network [19]. 

3.2 Data processing 

This analysis focuses on the search for (anti)muons, produced by muon-(anti) neutrino charged 
current interactions. Throughout the rest of this paper, muon-neutrinos denotes both muon- 
neutrinos and muon-antineutrinos. A data filter algorithm that is sensitive to muons from any 
direction has been applied during data taking. This algorithm first selects photons detected 
within 20 ns by separate optical modules in the same triplet. It then identifies a muon by 
requiring that for at least 5 triplets the relative arrival times of these photons are compatible 
with the signal expected from a muon traversing the detector. The muon purity for this data 
filter algorithm is better than 90%. Its average event rate was 1.0 Hz in the period considered 
and is mainly due to the background from atmospheric muons. 

These events were reconstructed offline to determine the muon trajectory, using a multi- 
stage fitting procedure. The reconstruction code follows the algorithm described in ref. [20]. 
Minor modifications were made to improve agreement between data and Monte Carlo. The 
most important modification compared to ref. [20] is that the amplitude information of the 
detected photons (hereafter referred to as hits) is discarded. As a result, the initial selection 
of signal hits is purely based on coincidences and the causality criterion. The algorithm 
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GRB 


GCN Circular 


GCN Report 


GRB 


GCN Circular 


GCN Report 


GRB070207 


6089 




GRB070724 




74.2 


GRB070209 


6101 


32.2 


GRB070729 




77.1 


GRB070227 




37.1 


GRB070805 


6711 




GRB070311 


6189 




GRB070808 


6724 




GRB070326 


6653 




GRB070810B 




81.1 


GRB070328 




42.3 


GRB070821 


6766 




GRB070330 




43.2 


GRB070824 


6768 




GRB070419 




47.1 


GRB070911 


6810 


82.2 


GRB070419B 




49.1 


GRB070913 




88.1 


GRB070420 




48.1 


GRB070923 


6821 




GRB070429B 


7140 


51.1 


GRB071001 




86.1 


GRB070508 




54.2 


GRB071003 




87.2 


GRB070509 




55.1 


GRB071008 




91.1 


GRB070517 




56.2 


GRB071010 


6864 


89.2 


GRB070611 




63.3 


GRB071010B 


6888 


92.1 


GRB070612 


6556 


64.1 


GRB071017 


6927 




GRB070612B 




65.1 


GRB071028B 




105.1 


GRB070615 


6537 




GRB071031 




99.1 


GRB070707 


6615 




GRB071112B 




103.1 


GRB070721 




72.2 


GRB071117 




106.2 



Table 1. List of the GRBs that are included in the analysis. The column "GRB" indicates the name 
of the GRB. The columns "GCN Circular" and "GCN Report" indicate the numbers of the GCN 
Circulars archive [15] and GCN Reports [16] respectively, from which the measured GRB data were 
taken. 



starts with a linear prefit which is used as a starting track for the subsequent stages. In 
addition, eight different starting tracks are generated by rotating and translating the result 
of the prefit. This is to increase the probability to find the global maximum of the likelihood 
function. The final stage of the fitting procedure consists of a maximum likelihood fit of the 
measured photon arrival times. A quality parameter for the fit, indicated by A, is determined 
based on the final value of the likelihood function. The A parameter is quantified by 

A = i° g(L) . + 0- 1 >< (^comp " 1), (3.1) 

which incorporates the maximum value of the likelihood, L, and the number of degrees of 
freedom of the fit, i.e. the number of hits, N^g, used in the fit minus the number of fit 
parameters; N comp is the number of times the repeated initial steps of the reconstruction 
converged to the same result. In general, N comp = 1 for badly reconstructed events while it 
can be as large as nine for well reconstructed events. The coefficient 0.1 in eq. 3.1 was chosen 
to maximise the separation in A between simulated signal and misreconstructed downgoing 
muons. The same reconstruction algorithm was used for the analysis described in ref. [21] 
which contains a brief description of the algorithm. For a more detailed description of the 
reconstruction algorithm see ref. [20]. 

The fit also provides an estimate of the uncertainties on the track parameters. These are 
used to select events with a well defined direction. The distribution of the estimated angular 
uncertainty on the direction of the reconstructed muon track, indicated by /3, is shown in 
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Figure 2. Left: Distribution of the estimated angular uncertainty (3 on the direction of the re- 
constructed muon track for upward reconstructed simulated atmospheric muons (dotted line), atmo- 
spheric muon-neutrinos x 1300 (dashed line) and signal events with an E~ 2 spectrum (dash-dotted 
line) normalised to the atmospheric muon-neutrinos, compared to upward reconstructed events in 
the data covering about 7 days and containing the 40 selected GRBs (data points), without fur- 
ther cuts. Right: Cumulative A distribution for upward reconstructed simulated atmospheric muons 
(dotted line), atmospheric muon-neutrinos (dashed line) and signal events with an E~ 2 spectrum 
(dash-dotted line) normalised to the atmospheric muon-neutrinos, compared to upward reconstructed 
events in the data covering about 7 days and containing the 40 selected GRBs (solid line), for events 
with (3 < 1°. The vertical lines indicate the analysis cuts. 

figure 2 (left). Since all selected GRBs occurred below the ANTARES horizon, only events 
with an upward reconstructed direction are considered. 

Due to the changing detector conditions during the period considered, and the varying 
environmental conditions, the statistics of the uncorrelated data that are equivalent to the 
conditions during the GRBs is too limited to estimate the background. Hence the background 
estimate has been made based on simulations. Figure 2 includes the expected distributions 
from atmospheric muons, atmospheric muon-neutrinos as well as an assumed signal with an 
E~ 2 spectrum, obtained from simulations. The neutrino signal is generated with software 
packages [22] that simulate the neutrino interaction as well as the production and propaga- 
tion of charged particles. The simulation uses the model for the atmospheric muon- neutrino 
flux from ref . [23] . The atmospheric muon contribution is simulated with the MUPAGE pack- 
age [24], which is based on a full Monte Carlo simulation of primary cosmic ray interactions 
and shower propagation in the atmosphere and reproduces the MACRO data [25, 26]. The 
simulated atmospheric muon contribution has an equivalent live time of one month. 

In the simulations, the stochastic energy loss of the muons, the production and prop- 
agation of the Cherenkov photons, the response of the photo-multiplier tubes to Cherenkov 
light and the simulation of the detector electronics are all included. The simulated photon 
signals are processed with the same data filter and reconstruction algorithms as the data. To 
obtain a realistic simulation of the varying environmental conditions due to bioluminescence, 
the measured optical background in coincidence with the prompt emission of each of the 
selected GRBs is taken from the data and added to the simulated events. 




Figure 3. Left: the angle averaged neutrino effective area of the detector consisting of 5 detection 
lines, and the effective areas for different neutrino zenith angle bins (cos (6) = 1 corresponds to ver- 
tically upward). The selection cuts are included. Right: the cumulative distribution of the angle 
between the reconstructed muon direction and the true neutrino direction for the detector consisting 
of 5 detection lines. The distribution is shown for selected events with an E~ 2 spectrum. 



Selection cuts on the A and /3 parameters are set to achieve a reliable rejection of misre- 
constructed atmospheric muons, while keeping a high signal efficiency. Neutrino candidates 
are required to be reconstructed as upgoing muons, with an estimated angular uncertainty 
/3 < 1°, see figure 2 (left). This cut removes 72% of the misreconstructed atmospheric muons. 
Figure 2 (right) shows the cumulative A distribution for upgoing reconstructed events, where 
the cut of (3 < 1° has been applied. Neutrino candidates are, in addition, required to have 
a quality value of A > —5.5. The cuts on the A and (3 parameters are chosen such that the 
background rate is reduced to a level below 2 x 10 -5 Hzsr -1 , while the signal efficiency is 
about 60% for the models considered in this analysis. The strong constraints from the time 
and direction coincidence with a GRB allow for a looser cut on the quality parameter A than 
that applied for the ANTARES point source search [27]. 

In both plots of figure 2, the contribution from atmospheric muons was scaled by 0.85 
to reproduce the data. This scaling factor is well within the uncertainty on the flux normal- 
isation [28]. 

The neutrino effective area of the detector consisting of 5 detection lines is shown in 
figure 3 (left) as a function of the neutrino energy. The presented effective areas include 
the selection cuts and are the average of the effective areas for muon-neutrinos and muon- 
antineutrinos. The angle averaged effective area is shown as well as the effective areas for 
different neutrino zenith angle bins. For vertically upward-going neutrinos the effective area 
is suppressed at high energies due to the absorption in the Earth. 

The angular resolution of the detector is determined by the angular uncertainty on 
the reconstructed muon direction (/3) and the neutrino scattering angle. The distribution 
of the angle between the reconstructed muon direction and the true neutrino direction for 
the detector consisting of 5 detection lines was evaluated with simulations. The cumulative 
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Figure 4. Top: the estimated background rate per unit solid angle as a function of the reconstructed 
muon direction in detector coordinates (left: zenith angle cosine 6, right: azimuthal angle <f> in radians; 
cos(6) = l corresponds to vertically upward). Bottom: the location of each of the 40 selected GRBs. 



distribution of this angle is shown in figure 3 (right) for events with a neutrino spectrum 
proportional to E~ 2 that passed the selection cuts. The median of this angular resolution 
was estimated to be 0.5 ± 0.1 degrees. 



4 Data analysis 

The remaining background due to atmospheric muon-neutrinos (89%) and misreconstructed 
atmospheric muons (11%), after the cuts on f3 and A, was estimated from simulations. The 
measured optical background from the data in coincidence with the prompt emission of 
each of the selected GRBs was added to the simulated events. The background rate per 
unit solid angle as a function of the reconstructed muon direction in detector coordinates 
is shown in figure 4 (top plots). The contribution from atmospheric muons was scaled by 
0.85 to reproduce the data. The background rate as a function of the reconstructed zenith 
angle, 6, is less isotropic compared to the background rate as a function of the reconstructed 
azimuthal angle, 4>. However, since the selected GRBs are distributed rather isotropically 
(bottom plots of figure 4), the average background rate may be used for each GRB. The 
expected background rate per unit solid angle is estimated by taking a weighted average of 
the background rate in each solid angle bin, where the weighting accounts for the relative 
duration of GRBs in that solid angle bin. This results in an estimated background rate of 
1.54 x 10" 5 Hzsr" 1 . 



- 8 - 



A neutrino candidate is considered to come from the GRB when the reconstructed 
muon track points back to the GRB within 2°. Prom simulations it was estimated that for 
about 85% of the neutrino candidates, the muon is reconstructed within 2° from the neutrino 
direction. This high signal efficiency is the result of the selection of events with a good angular 
resolution (the cut on the /3 parameter). The uncertainty on the source positions of the GRBs 
considered (table 1) is smaller than 0.07°. The total background, fibkg, due to atmospheric 
muon-neutrinos and misreconstructed atmospheric muons during the prompt emission of the 
40 GRBs can be expressed as 

dR 

Vbkg = ^ X X T ( 4 ' 1 ) 

where ^ is the background rate per unit solid angle, is the solid angle of the 2° search 
cone and T is the total prompt emission duration of the 40 GRBs. For the prompt emission 
duration of the GRB, the so called Tgo time interval is taken [29], extended by 5% before 
the start, and 5% after the end of this time interval. T amounts to 2114 seconds. The total 
estimated background amounts to 1.24 x 10 -4 events. 

5 Results 

This analysis focuses only on the detection of muon-neutrinos in coincidence with the prompt 
emission phase of the GRB. The search for neutrinos in correlation with GRBs was done in 
a stacking analysis, in which all data in coincidence with the prompt emission of the GRBs 
were accumulated. 

Selected neutrino candidates are considered to be correlated with a GRB when their 
detection time is in coincidence with the prompt emission of the GRB, assuming that neu- 
trinos travel at the speed of light. After unblinding the data, no neutrino candidates were 
found in correlation with the selected GRBs. 

5.1 GRB neutrino spectra 

A general neutrino spectrum of E~ 2 is assumed for the neutrino emission from GRBs, where 
E u is the neutrino energy. In addition, three other energy spectra for neutrino emission from 
GRBs have been considered: the energy spectrum proposed by Guetta et al. [30], the general 
Waxman and Bahcall energy spectrum [31], and the energy spectrum proposed by Ahlers et 
al. [32]. 

The energy spectrum according to Guetta et al. is calculated for each GRB individually 
using the data from the instruments on the satellite that detected the GRB, taken from the 
references given in table 1. The method to calculate a neutrino fluence from an individual 
GRB is provided in ref. [30]. The prediction therein referred to as "Model 2" is used. In 
case not all required parameters were measured, nor default values are provided in ref. [30], 
the values as listed in table 2 were used. Within the list of GRBs no single GRB yields 
a detectable signal. The total estimated muon-neutrino fluence for the 40 selected GRBs 
according to Guetta et al. is the sum of the calculated individual muon-neutrino fluences, 
and is shown in figure 5 (left). The total number of expected events from the 40 selected 
GRBs for the estimated muon-neutrino fluence with the energy spectrum according to Guetta 
et al. is 1.7 x 10 -3 . The estimated muon-neutrino fluence represents the fluence at Earth and 
includes the effect of neutrino oscillations. 

Waxman and Bahcall [31] provide an approximate estimate of the neutrino energy spec- 
trum, which is assumed to be the same for each GRB. The expected number of events for 40 
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Parameter 


Value 


z 


0.25 or 2.8 


e e 


0.1 


Q>> 


0.1 


J> 

C 7 


0.4 MeV 


a ~ 


1 




2 




10 51 ergss -1 




6 x 10 -6 ergs cm -2 



Table 2. Assumed parameter values used for calculating the neutrino spectrum for an individual 
GRB using the method described in rcf. [30] (see text). The symbols correspond to those used in 
rcf [30]. The rcdshift z is set to 0.25 for GRBs with a prompt emission duration of less than 2 s [33], 
and set to 2.8 for other GRBs [34]; the values for the fraction of the internal energy in electrons 
(e e ) and the magnetic held (cb) are taken from [35]; the default values for the break energy in the 
j-ray spectrum e^, the spectral indices of the j-ray spectrum before (ot^) and after (f3 7 ) er,, the 7-ray 
luminosity L 1 , and the 7-ray fiuence Ky are the mean values from the Swift catalogue [36, 37]. 
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Figure 5. Left: the 90% CL upper limit on the muon-neutrino fiuence F v from the 40 selected 
GRBs for the energy spectrum according to Guetta et al. (solid line) and the corresponding total 
estimated muon-neutrino fiuence including oscillations (dashed line). Right: the 90% CL upper limit 
on the diffuse muon-neutrino flux for the E~ 2 energy spectrum (thick solid line), the Waxman 
and Bahcall energy spectrum (thin black line), and the energy spectrum according to Ahlers et al. 
(thin grey line). The black dash-dotted line is the 90% CL upper limit of IceCube [2] assuming the 
Waxman and Bahcall energy spectrum. Also shown are the estimated diffuse muon-neutrino fluxes 
including oscillations assuming the Waxman and Bahcall energy spectrum (black dashed line) and 
that according to Ahlers et al. (grey dashed line). 



GRBs assuming the Waxmann and Bahcall energy spectrum is 7.0 x 10~ 3 . The estimated 
GRB muon-neutrino intensity assuming a Waxman and Bahcall energy spectrum is presented 
as a diffuse flux in figure 5 (right) where a total of 10 3 GRBs are expected per year. 1 The 



1 In ref. [31] the evolution correction is assumed to be 1. 



right plot in figure 5 also includes an estimated diffuse muon-neutrino flux according to 
Aiders et al. [32] where typical values for the parameters of the GRB environment were used 
as provided in ref. [32]. The estimated diffuse muon-neutrino fluxes for both energy spectra 
represent the fluxes at Earth and include the effect of neutrino oscillations. 

5.2 Upper limits on the diffuse neutrino flux and the neutrino fluence 

Since no events were found in correlation with the prompt photon emission of the 40 selected 
GRBs, upper limits have been placed on the intensity of the diffuse muon-neutrino flux 
and the muon-neutrino fluence at Earth for the different models. The 90% confidence level 
limits were set using the Feldman-Cousins prescription [38], and are shown in figure 5. The 
same systematic uncertainties as described in ref. [27] have been considered. These include 
the effect of reduced optical module efficiencies and the effects which have a net result of 
degrading the time resolution, such as possible mis-alignments of the detector, inaccuracies 
in the simulation of light propagation in the water or in the transit time distribution of 
the photo-multiplier tubes. The impact of the systematic uncertainties were evaluated by 
including these effects in the simulation, described in ref. [27]. This results in a degradation 
of the limits of less than 10%. 

For the E~ 2 energy spectrum, 90% of the signal is expected in the energy range 
5.2 TeV< E v <1.4 PeV. The upper limit on the diffuse muon-neutrino flux for the E~ 2 
energy spectrum is 2.7 x 10~ 7 (E v /Ge\)~ 2 GeV" 1 cm" 2 s" 1 sr" 1 . For the energy spec- 
trum according to Waxman and Bahcall 90% of the signal is expected in the energy range 
41 TeV< E v <2.1 PeV and a 90% CL upper limit of about 294 times the predicted flux 
as shown in figure 5 (right) was set. For the energy spectrum according to Ahlers et al. 
90% of the signal is expected in the energy range 33 TeV < E u < 2.0 PeV and a 90% CL 
upper limit of about 11 times the predicted flux as shown in figure 5 (right) was set. For 
the energy spectrum according to Guetta et al. 90% of the signal is expected in the energy 
range 86 TeV < E v < 6.0 PeV and a 90% CL upper limit of about 1467 times the predicted 
fluence as shown in figure 5 (left) was set. 

6 Conclusions 

A search for muon-neutrinos in correlation with the prompt emission of gamma-ray bursts 
using the data taken with the ANTARES detector during the first year of operation has been 
presented. During the period considered, the detector was less than half its final size. The 40 
GRBs that were examined for neutrino emission were selected from the GRB observations by 
satellite instruments. No correlations between neutrinos and the selected GRBs have been 
found. Upper limits have been obtained on the fluence of muon-neutrinos from the 40 GRBs 
and on the diffuse muon-neutrino flux for different models. 

A neutrino telescope in the Mediterranean Sea is well suited to detect high-energy neu- 
trinos in correlation with GRBs spread over a wide region of the sky, including in particular 
the Southern hemisphere. The low background is the result of the short duration of GRBs 
and the excellent angular resolution. The ANTARES detector was completed mid 2008 and 
is 2.5 times bigger than the detector configuration considered in this analysis. Since the com- 
pletion of the detector, on average 250 GRBs per year have been detected in the Southern 
hemisphere and at least five more years of data taking are foreseen. With this, a large sample 
of GRBs is available for further analysis, complementary to the IceCube field of view and en- 
ergies. The next-generation neutrino telescope KM3NeT [9], to be built in the Mediterranean 
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Sea, will surpass the ANTARES sensitivity by two orders of magnitude. Once operational it 
will probe GRB models with unprecedented sensitivity, and will push the boundaries towards 
new discoveries. 
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